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Abstract. By using the corrections to the asymptotic scaling forms of 
the fields of the O(N) Gross Neveu model to solve the dressed skeleton 
Schwinger Dyson equations, we deduce the critical exponent corresponding 
to the /^-function of the model at 0(1/N 2 ). 
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1 Introduction. 



One of the persistent problems of quantum field theories is a lack of total 
knowledge of the renormalization group functions, such as the /3-function, 
which are important for having a precise picture of the quantum structure 
of a theory. For certain models the functions can be computed to three or 
four orders as a power series in a perturbative coupling constant which is 
assumed to be small. However, in explicit calculations one has to work ap- 
preciably harder to gain new information which is partly due to the increased 
number of Feynman graphs one has to analyse within some renormalization 
scheme such as MS, even in theories with the simplest of interactions. It 
is therefore important to develop different techniques to give an alternative 
picture of the perturbative series. One method which achieves this is the 
large N expansion for those theories where one has an iV-tuplet of funda- 
mental fields. Then 1/N is a small quantity for N large and this can be used 
as an alternative expansion parameter. Whilst a conventional leading order 
analysis is relatively straightforward to carry out for most theories it turns 
out that it is not useful for going to subsequent orders in 1/N. To obviate 
these difficulties methods were developed for the O(N) a model which were 
successful in solving that model at 0(1/N 2 ), [||, §]. In particular the method 
uses a different approach to the conventional renormalization of the large iV 
expansion in that one solves the field theory precisely at its (i-dimensional 
critical point, ie the non-trivial zero of the /^-function, by solving for the 
critical exponents. Performing the analysis at the fixed point of the renor- 
malization group means that there are several simplifying features. First, 
the theory is finite. Second, since (5{g) is zero the theory has a conformal 
symmetry and the fields are massless. This has two consequences, one of 
which is that the critical exponents of various Green's functions can be de- 
termined order by order in 1/N in arbitrary dimensions, Jl], gj. Second, 
the masslessness of the fields simplifies the Feynman integrals which occur 
and allows one to compute the graphs in arbitrary dimensions, which are 
otherwise intractable in the conventional (massive) large N renormalization. 
By solving for the exponents in this fashion one can then relate the results 
through an analysis of the renormalization group equation at criticality to 
the critical renormalization group functions. (See, for example, Q.) Hence, 
one gains, albeit by a seemingly indirect method, information on the pertur- 
bation series of the theory to all orders in the coupling at the order in 1/N 
one is interested in via the e-expansion of the exponent. Clearly, this has 
important implications for gaining a new insight into the renormalization 



2 



group functions at large orders of the coupling as well as allowing one to 
check the series with explicit calculations at low orders. Furthermore, since 
one calculates in arbitrary dimensions a three dimensional result will always 
be determined simultaneously. 

Since the earlier work of [jl], the method has been extended to models 
with fermions, j|, ||, supersymmetry || 0] and theories with gauge fields 
[8-11]. In this paper we present the detailed evaluation of the /^-function ex- 
ponent for the four-fermi theory or the O(N) Gross Neveu model, | |12| |. The 
motivation for such a calculation is that first of all knowledge of 2 A = — (3'(g c ) 
at 0(1/ N 2 ) will mean that the field theory will be solved completely to this 
order. In Q the exponent rj, which is the fermion anomalous dimension, 
was calculated at 0(1/N 2 ) and more recently the vertex or vp-ip anomalous 
dimension was also computed to the same order, Q. Together with A one 
can deduce the remaining thermodynamic exponents for this model through 
the hyperscaling laws discussed in [13|. Secondly, the computation of A- 
type exponents for theories with fermions as their fundamental field is not 
as straightforward as the case where one deals with purely bosonic fields. As 
was noted in [11] there is a subtle reordering of the graphs in the formalism 
and it is important to have a complete understanding of this feature if one is 
to apply similar methods to deduce results in physical gauge theories. Also 
to a lesser extent the methods which we had to develop to solve the current 
problem, which essentially is the evaluation of massless four loop Feynman 
diagrams will prove to be extremely useful in other contexts. Finally, we are 
interested in going well beyond the leading order in the three dimensional 
Gross Neveu model to compare estimates of various exponents from our an- 
alytic work with Monte Carlo simulations currently being carried out, [14]. 
The leading order results are not precise enough to be able to compare with 
the relatively low values of N which are being simulated and the 0(1/ N 2 ) 
reults therefore must be computed. 

The paper is organised as follows. In section 2, we introduce our nota- 
tion and review the leading order formalism used to compute the exponent 
A at 0(1/ N) which will serve as the foundation for the 0(1/N 2 ) corrections. 
This formal extension is discussed in section 3 where we derive finite consis- 
tency equations and explain the need to compute some three and four loop 
Feynman graphs. The explicit evaluation of these is discussed in sections 4 
and 5 whilst the 0(1/ N 2 ) corrections to a 2-loop integral which appears at 
0(1 /N) are derived in section 6. We conclude our calculation in section 7 
by giving an arbitrary dimensional expression for A at 0(1/N 2 ) and discuss 
the numerical predictions deduced from it for the three dimensional model. 
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2 Preliminaries. 



The theory we consider involves self interacting fermions if) 1 , 1 < i < N, 
where 1/N will be the expansion parameter for N large. One can formu- 
late the lagrangian either by using the explicit four point interaction or by 
introducing a bosonic auxiliary field, a, which is the version we use. The 



quantum theory of both are equivalent. Thus we take, [12] 



L = IfiW + \aff - ^ (2.1) 

where g is the perturbative coupling constant which is dimensionless in two 
dimensions. The aim will be to calculate the 0(1/N 2 ) corrections to the 
/3-function and we note that the three loop structure of this has already 
been calculated perturbatively in dimensional regularization using the MS 



scheme, ]15|, 16[ [17|, 18|, as 



(3(g) = (d-2)g-(N-2)g 2 + (N-2)g 3 + \(N-2)(N-7)g 4 (2.2) 

where the coupling constant in (2.2) is related to that of (2.1) by a factor 
2tt which we omit here since it will play a totally passive role in the rest 
of the discussion. It is important to note that in carrying out perturbative 
calculations with dimensional regularization that (2.2) is what one deter- 
mines as the /3-function in ci-dimensions prior to setting d = 2 to obtain the 
renormalization group functions in the original dimension. There, of course, 
the theory is asymptotically free. However, the (i-dimensional /3-function 
(2.2) can be viewed from a different point of view for the large N critical 
point analysis of the present work. For instance, when d > 2 which is the 
case we will deal with for the rest of the paper, there exists a non-trivial 
zero of the /3-function at a value g c given by 

9c ~ t — - — r ( 2 -3) 
y (N-2) 17 

at leading order in large N where the corrections are 0(e 2 ) and 0(1/ (N— 2) 2 ) 
and d = 2 + e. This corresponds to a phase transition which is apparent in 
the explicit three dimensional work of ||i~9|, pofl. Indeed similar approaches 



were examined in the work of [21] for the 0(N) a model. When one is in the 
neighbourhood of a phase transition, it is well known that physical quan- 
tities possess certain power law behaviour. For physical systems the power 
or critical exponent fundamental to the power law totally characterizes the 
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properties of the system. In continuum field theory, in the neighbourhood of 
a phase transition, the Green's functions also exhibit a power law structure 
where the critical exponent, by the universality principle, has certain prop- 
erties, J3j. For instance, it depends only on the spacetime dimension and 
any internal parameters of the underlying field theory. More importantly, 
though, for our purposes, one can solve the renormalization group equation 
at criticality and relate several exponents to the fundamental functions of 
the renormalization group equation, H, which like (2.2) are ordinarily cal- 
culated order by order in perturbation theory. In the alternative critical 
point approach one can compute the exponents at several orders in 1/N 
which then gives independent information on that critical renormalization 
group function. Since the location g c is known at leading order as a func- 
tion of e and 1/N, one can undo the relations between exponent and critical 
renormalization group function to deduce the coefficients appearing in the 
perturbative series. Clearly this is a powerful alternative method of com- 
puting, say, /3-functions. 



For this paper, we extend the earlier work of [11] which was based on the 
pioneering techniques developed for the bosonic a model on S N , [|], [2|. The 
physical ideas behind the method are relatively simple. In the neighbour- 
hood of g c the model is conformally symmetric and therefore the Green's 
functions scale. To analyse the critical theory one postulates the most gen- 
eral structure the Green's functions can take which is consistent with Lorentz 
and conformal symmetry, Q . The critical exponents of these scaling forms 
involves two pieces. One is related in the case of a propagator to the canoni- 
cal dimension of the field as defined by the fact that the classical action with 
lagrangian (2.1) is a dimensionless object. Since quantum fluctuations will 
always alter the canonical dimension a non-zero anomalous dimension is ap- 
pended to the canonical dimensional and it carries the information relevant 
for the renormalization group functions. 

To be more concrete and to fix notation, for (2.1) the scaling forms of 
the propagators in coordinate space as x — > are, 



Ad B 

~ jjy > < x ) ~ J^y ( 2 - 4 ) 



where 

a = \x + \r\ , p = 1 - 7] - x (2.5) 

and rj is the fermion anomalous dimension, x is the vertex anomalous di- 
mension and d = 2fi is the spacetime dimension. Both r] and x have been 
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calculated at 0(1/N 2 ) within the self consistency approach, [||, ||, and are 

2(^-i) 2 r(2^-i) 



m 



Vi 



2^ -If 



A* 



r(2-/i)ro* + i)r2(/i) 

+ 3/i + 4(/i - 1) + 2(^x - 1)(2 M - 1)^(m) 
Wl 



Xi 



(m-i) 



(2.6) 

(2.7) 
(2.8) 



X'2 



ml 



01-1)3 

(2^-l)(^ 2 



[3/x(// - 1)0(^) + (2// - l)^(^) 



M-l) 



0*-i) 



(2.9) 



where = tfj(2fi - 1) - ^(1) + if>(2 - m) - ^0«), @(a0 = tf(p) ~ ^(1) 
and ^(A 4 ) is the logarithmic derivative of the T-function. The expression 
(2.6) was first derived in ||| and later in [23-25] where xi was & l so given 
in pSj ]. Each expression (2.6)-(2.9) agrees with the respective three loop 
perturbative results of the corresponding renormalization group functions 
in d = 2 + e dimensions which were given in [16]. The quantities A and -B 
in (2.4) are the amplitudes of ip and a respectively and are independent of 
x. 

The method to deduce (2.6) and (2.7) is to take the ansatze (2.4) and 
(2.5) and substitute them into the skeleton Dyson equations with dressed 
propagators of the 2-point function, EL ffl, which are valid for all values of 
the coupling including the critical coupling. One subsequently obtains a 
set of self consistent equations which represent the critical Dyson equations. 
Their solution fixes r\i which is the only unknown at the ith order where 
V = Si^i Vi/N 1 - Further the vertex anomalous dimension x is determined 
by considering the scaling behaviour of the aipip vertex also in the critical 
region using a method developed in || which extended the earlier work of 
@ to 0(1/N 2 ). 

To determine the corrections to the /3-function one follows the analogous 
procedure used in ]2| and developed for models with fermion fields in [Q, p"l| ]. 
If we set 2A = —(3'{g c ) then the critical slope of the /3-function can be 
computed by considering the corrections to the asymptotic scaling, M , ie 



V>(x) 



A 



2\A 



1 + A'{x 2 ) 



6 



a(x) 



B 



1 + B'l 



„2\A 



(2.10) 



where A' and B' are new amplitudes and A = ^ — 1 + Yli^i from 
(2.2). The idea then is to compute A at 0{\/N 2 ) in arbitrary dimensions. 
Once obtained we can use the relation between A and (i'{g c ) to deduce (3{g) 
as a power series in g at the same approximation in large N, since knowledge 
of Ai allows us to determine the value of g c to undo the relation. 

We close the section by reviewing the method of [||, |ll]] to deduce Ai. 
As indicated we use the skeleton Dyson equations which are illustrated in 
fig. 1. To deduce r\\ the equations were truncated by including only the one 
loop graphs of fig. 1. However, it turns out, as we will recall below, that 
for Ai one has to consider the additional two loop graph of the a equation. 
The quantities tp -1 and cr _1 are the respective two point functions and 
their asymptotic scaling forms have been deduced from (2.4) by inverting in 
momentum space using the Fourier transform 



1 



a{a) 



where a(a) 



V^Or) 
a-\x) 



(x 2 ) a 2 2a irv J k {k 2 y- a 
a)/T(a). Thus as x — » 0, M, 
r{a — l)^f 



A(x 2 ) 2 ^- a+1 



1 



1- A's(a-l)(x 2 ) x 
-B'q((3)(x 2 ) x 



where 



a(P - (J,) 



r(a) 



vr 2 ^a(/3) 
a(P-H + \)a(/3 - A) 
a{[5 - /i)a(/3) 



ap(a) 
, s(a) 



(2.11) 

(2.12) 
(2.13) 

(2.14) 



(a — fi + A) (a — A) 



To represent the graphs of fig. 1 one merely substitutes (2.10), (2.12) and 
(2.13) for the lines of each of the graphs to obtain 








r(a - 1)[1 - A's(a - l)(x 2 ) x ] + z[l + (A' + B')(x 2 ) x ] (2.15) 



p{0) 



(x 2 ) 2 ^-! 3 
Nz 2 



-[l-B' q (P)(x 2 ) x } + 



Nz 



(x 



2\2a- 



T [1 + 2A'(: 



\n 1 + (TL 1A A , + U 1B B')(x 



2\X] 



(2.16) 
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where we have not cancelled the powers of x 2 in (2.10) and the quantities 
IIi, IIi^ and Uib are the values of the two loop integral in the respective 
cases when there are no (x 2 ) x contributions, when (x 2 ) A is included on a ifi 
line and when it is included on the a field. We have also set z = A 2 B. As in 
[^, |[] the terms of (2.15) and (2.16) involving powers of (x 2 ) x decouple from 
those which do not to leave two sets of consistency equations. One set yields 
r\\ whilst the second determine Ai. To achieve this one forms a 2 x 2 matrix 
which has A! and B' as the basis vectors and sets its determinant to zero 
to have a consistent solution. It is the subtlety of taking this determinant 
which necessitates the inclusion of YI\b, H], in (2.16) as discussed in [11]. 



Basically when one substitutes the leading order values for a and f3 into the 
basic functions (2.14) one finds 

s(a-l) = 0(N) , r(a-l) = O ( -M , q((3) = o(±-\ (2.17) 



NJ ' * v ' \N. 

Thus analysing the leading order N dependence of each of the elements of 
the matrix one finds that the contribution from the terms in a^ 1 involving 
B' are of the same order as the (finite) two loop graph Uib- Thus it cannot 
be neglected and we note that explicit evaluation gave 

UlB = M)¥p (2 ' 18) 

Thus substituting into the equation 

/ — r(a — l)s(a — 1) z \ . 

det( ^ _mm.^ 1B ) - < 2A9 > 

one deduces 

Xi = - (2/i - l)r/i (2.20) 

as was recorded in Q . This completes our review of the previous work in this 
area and lays the foundation for the subsequent higher order calculations. 



3 Master equation. 

In this section, we derive the formal master equation whose solution will 
yield A2- As already indicated in the previous section this involves trun- 
cating the Dyson equations at the next order and including the appropriate 
corrections. For the moment we concentrate on the equation for tp as it has 
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a simpler structure compared to (2.16). The additional 0(1/N 2 ) correction 
we consider is illustrated in fig. 2 and we denote it by E. Including it in 
(2.15) we have 

° = $2~li [1 ~ M« - 1)(* 2 ) A ] + pyg^ZA [1 + ^ + ^)(* 2 ) A ] 



-2a-l-2A P + + ^B)(X 2 ) A ] (3.1) 



where the subscript on the corrections E^ and E# correspond to the inser- 
tion of (x 2 ) x on either the ip or a lines of the graph of fig. 2 and E, E^ 
and T,b are the values of the respective integrals. There are two graphs 
making up Eg due to the presence of two a lines and each give the same 
contribution. For E^, there are three graphs, one of which gives a different 
value from the other two where the insertion is on a line adjacent to the 
external vertex. In (3.1) we have included the additional quantities A and 
m. The graph E arises in Q in the determination of r/2 and it is in fact 
infinite which can be seen by the explicit computation using the uniqueness 
method developed first in [27] and later in g [2§]. Consequently one has 



to introduce a regularization by shifting the exponent of the a field by an 
infinitesimal quantity A, ie (3 — > (3 — A. To remove the infinities from E, 
Tia and E# one uses the counterterm available from the leading order one 
loop graph. Thus formally setting 

S = f + E' , E A , B = Zjj* + X> A>B (3.2) 
in (3.1) and expanding 

"* = 1 + H + °(^) (3 ' 3) 

the divergent terms of (3.1) are set to zero minimally to obtain a finite 
consistency equation ie 

ziK A ziK B 
-i = - — = - — (3-4) 

which implies K A = Kb and this will provide a check on the explicit cal- 
culation described later. In order to proceed to the critical region one must 
exclude the lnx 2 style terms which remain which is achieved by exploiting 
the freedom in the definition of the vertex anomalous dimension by setting 

Xi = - z 1 K A = - ziK B (3.5) 
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Agreement with (2.8) will be another check. This will leave a finite set of 
equations which are valid as x — > which again decouples into one which is 
relevant for 772 and the other for A2 ie 

= r(a-l) + z + z 2 Yl (3.6) 

from which Z2 can be derived and 

= A'{z-r(a-l)s(a-l) + z 2 Z' A } + B'[z + z 2 Z' B ] (3.7) 

However, by analysing the iV-dependence of each term of the A' coefficient 
of (3.7) the correction is 0(l/N 2 ) with respect to r(a — l)s(a — 1) and 
therefore does not need to be computed explicitly since it will contribute to 
A3 and not A2. 

We now turn to the a equation. In the same way we had to consider the 
higher order two loop graph of fig. 1 to deduce Ai , we now have to include 
the analogous set of graphs for the next order to determine A2. As we are 
using graphs with dressed propagators it turns out there are only five graphs 
which arise. These are illustrated in figs 3 and 4 and we have given each 
a label. The subscript B indicates that we need only consider the graphs 
where there is an (x 2 ) A insertion on the a line. Again the insertions on the ifi 
lines will be relevant for A3. We have grouped the graphs which are divergent 
and therefore require regularization by A. The origin of the infinity is the 
same as the vertex infinity which occurs in E. Indeed in fig. 3 each graph 
corresponds to the usual vertex correction of the two loop graph IIi. For 
the A consistency equation of the a Dyson equation there is an insertion of 
(x 2 ) A in one of the a lines of the graphs. For the case II3B, for example, the 
insertion in one line removes the infinity from that vertex since the presence 
of the exponent A on the a lines moves the overall exponent of that line 
from the value which gives infinity. Thus the graph has a simple pole in A 
which is removed by the same vertex counterterm as (3.4). For H2B one of 
the insertions onacr line makes the graph finite and we call it Il2_B2 and no 
regularization is required. The other case, H2B1, is divergent but is again 
rendered finite by (3.4) in the consistency equation. By contrast the graphs 
of fig. 4 do not involve any divergent vertex subgraphs and when any one 
of the a lines has an (x 2 ) x insertion each graph is completely finite. For 
notational convenience we define H5B1 to be the graph with an insertion on 
the top a line and II552 to have an insertion on the central a line. Similarly, 
we denote the graph of IlgB with the bottom a line corrected by Hqbi and 
the other case by Il6£2- Therefore, there are five distinct finite massless 
Feynman graphs to evaluate. 
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We have dealt at length with the higher order graphs which need to 
be included in the corrections to (2.16). However, there are also 0(1/N 2 ) 
contributions coming from H\b- For, Ai, one considers this graph with a 
= fj,, (3 = 1 and A = \i — 1. These are only the leading order values of 
the exponents and since we are dealing with fields with non-zero anomalous 
dimensions which are 0(1/N) these give contributions to A2 when J1\b is 
expanded in powers of 1/N. Therefore, these ought not to be neglected in 
deriving the master equation for A2. 

Rather than reproduce the analogous renormalization of the a consis- 
tency equation, which proceeds along the same straightforward lines as dis- 
cussed earlier, we obtain the following finite equation which includes the 
corrections to (2.16) 



= 2zA' - B' 



N + ~~2 + ~2 Um 



(3.8) 



where 



n B2 = 2n 2B1 + 2n 2B2 + 2n 3B + 2114B 

- z 1 [2n 5B1 + n 5i j 2 + 2n 6m +4n 6i j 2 ] (3.9) 

and the prime is understood on the integrals which are divergent. We have 
preempted the explicit calculation of later sections by using the fact that 
n' lyl = in writing down (3.8) and ignoring the corrections where there is an 
insertion on the tp lines of the graphs of figs. 3 and 4 since they are relevant 
for A3. 

This completes the derivation of the formal consistency equations which 
yield A 2 . One again sets the determinant of the 2x2 matrix formed by A' 
and B' as the basis vectors in (3.7) and (3.8) to zero, and all that remains 
is the evaluation of H\b and IIb 2 . 



4 Computation of divergent graphs. 

In this section we discuss the computation of the divergent graphs H\a, Sia, 
n 2 _Bi and H^b- First, though we recall the basic tool we use for computing 
massless Feynman graphs, which is the uniqueness construction first used in 



and developed for large N work in |2| and other applications in [28]. The 
basic rule for a bosonic vertex which we require is illustrated in fig. 5 whilst 
that for a a'lpip type vertex is given in fig. 6, 0]. In each case the arbitrary 
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exponents CKj and are constrained to be their uniqueness value, J2i a i = 
2ji and = 2/x + 1, whence the integral over the internal coordinate 

space vertex can be completed and is given by the product of propagators 
on the right side represented by a triangle. The quantity 1/(01,02,03) is 
defined to be 71^ nf=i «(««)• 

It is easy to observe that for the atpip vertex of (2.1) we have 2a + (3 = 
2/i + 1 at leading order so that in principle the integration rule of fig. 6 can 
be used. However, if we recall that there is a non-zero regularization A this 
upsets the uniqueness condition. To proceed with the determination of the 
divergent graph one instead uses the method of subtractions of Pj. Since 
we need only the simple pole with respect to A and the finite part of each 
divergent graph one subtracts from the particular integral another integral 
which has the same divergence structure but which can be calculated for non- 
zero A. The difference of these two integrals is A-finite and therefore can 
be computed directly by uniqueness. To determine the two loop graphs IIi^ 
and £1,4 which occur for A2 we refer the interested reader to the elementary 
definition of the subtracted integrals given in Q since the treatment of the 
three loop graphs is new and will be detailed here. 

First, we consider the case H^b- The subtraction we used is given in fig. 
7 where the right vertex subgraph is divergent and the subtraction is given 
by removing the internal tp line to join to the external right vertex. This 
graph can be computed for non-zero A and for a, (3 and A given by their 
leading order values. After integrating two chains one is left with the two 
loop integral 

(fj,, p,, /I - A, fi, 2 - fi) (4.1) 

where the general definition is given in fig. 8. To compute (4.1) we first 
make the transformations /" and \ in the notation of @ which leaves the 
integral (1, fl, fl,l — A, /z — 1) which has been computed already in Q. Thus 
overall we have 



. 27T 2/i 

(fi,p,,H- A,/i,2-/i) 



(M-i) 2 r 2 ( M ) 



l-£(3( M -l)e + 



2\ vr ' (fi-1) 



(4.2) 



To determine the remaining finite part one uses fig. 6 and as an intermediate 
step a temporary regularization 5 has to be introduced to perform integra- 
tions in both graphs in different orders, as in @, |, §. Useful in obtaining 
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the correct answer is the result 



(/x, /}, /}, p,, 2 - y, - A) = 



2vr 2 ' 1 



1 3Qi - 1)A 



G + 



0* - 1) ; 



(4.3) 



computed in a similar fashion to (4.2). After a little algebra the sum of the 
finite piece and subtracted integral yield 



n 3 s 



27T 4 " 



(/z-i)3rV)A 



A(/i-l) 



30 + 



1 



0* - 1) 2 



(4.4) 



To determine n 2 Bi we have illustrated one of two possible subtractions in 
fig. 7. The procedure is the same as that for II^b and also makes use of 
(4.2) and (4.3). We obtain 



II2B1 = - 
Finally, we note 

Sis = 



2vr 4 ^ 



(M-ifr^A 



i-A(M-i)(se + ' ) 



( M -i)r 2 ( M )A 



8tt 2 ^ 



(M-i)r 2 (/x)A 



(4.5) 



(4.6) 



where the finite part is zero in both cases and (4.6) are consistent with our 
choice of xi m the renormalization of the previous section. 



5 Computation of finite integrals. 

The remaining higher order graphs are A-finite and therefore do not need to 
be regularized. Moreover, we need only compute them for the leading order 
values of a, (3 and A. In other models, the higher order graphs were computed 
by uniqueness and we used this technique extensively for the calculation 
though we had to employ some novel methods which deserve discussion. 
As the integrals we need to determine involve massless propagators one 
can introduce conformal changes of variables on the internal vertices on 
integration. For example, one conformal transformation is 
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from which it follows that 

r 2 



x 2 — > -= (5.2) 

x z 



which was used in 0] to compute two loop graphs. For an integral with 
fermions the analogous situation for the propagators is 

* — i ( 5 - 3 ) 

from which we deduce that for a fermion propagating from x to y where 
both are internal vertices 

This latter transformation provides the starting point for computing each 
integral as it allows us to carry out several integrations over the internal 
vertices immediately. 

For example, we consider the three loop graph of fig. 4. First, integrating 
one of the unique vertices and then making a conformal transformation on 
the subsequent integral, which requires (5.1)-(5.4), one ends up with the 
first graph of fig. 9. The fermion trace is taken over the endpoints of 
the propagator with exponent fi joining to the vertex with two bosonic 
propagators and the propagator with exponent 1. The bosonic triangle 
of this graph is unique and can be replaced by a unique vertex. After 
integrating several chains one is left with the two loop graph of fig. 9. The 
techniques to compute two loop graphs are elementary. One performs the 
fermion trace which yields a series of chains of integrals and another integral 
proportional to the basic integral ChT(l, 1) defined in |2| as ChT(a\, 02) = 
(fj, — 1, a.\, a 2 , H — 1, n — 1) and evaluated as 



ChT(ai,a2) 



vr 2 ^a(2|i - 2) 

r(//-i) 



a(a%)a(2 — a%) 
(1 - a 2 )(«i + a 2 - 2) 



+ 
+ 



a(o2)a(2 — 02) 
(1 - ai)(«i + a 2 - 2) 
a(«i + «2 — l)a(3 — Qi — CK2) 



(a x - l)(a 2 - 1) 



Thus, 



n 4B 



(^-i) 2 r^) 



39 + 



(M-l) 



(5.5) 



(5.6) 
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For the four loop graphs of fig. 4, the conformal transformations (5.1)- 
(5.4) are again the starting point. For instance, in Uqbi the location of 
the exponent (2 — (a) and the topological structure of the graph means that 
elementary integrations quickly result in an integral of the form of fig. 10, 
where there is a fermion trace over the two propagators with both exponents 
1 and another over the remaining four propagators. In fact the unique vertex 
present in fig. 10 means that one only has to compute one two loop integral. 
This is again achieved by taking the fermion trace and one finds 



vr 6 ^a(2/i 



(M-l) 5 ^) 
7r 6 ^a 2 (2/i - 2) 



1 5 (2fi- 



Gu-1) 2(^-1)2 Ou-1) 



(5.7) 



The remaining three integrals, however, turned out to require a significant 
amount of effort. We consider U5B2 first. After several integrations one is 
left with the first graph of fig. 11 where again there are two fermion traces, 
one of which is over the two propagators with exponents 1. Taking this trace 
explicitly yields three graphs. Two of these are equivalent after several chain 
integrations and are proportional to the basic integral tr G(2 — fj,, 1) which 
was defined and evaluated in J/]]. The remaining integral is equal to the 
second graph of fig. 11, after performing the transformations <— on the 
right external vertex. Again taking the fermion trace yields two graphs 
which are equivalent and elementary to compute as they are proportional 
to ChT(2 — /u, 1) and a purely bosonic integral which is the third graph of 
fig. 11. It is completely finite. However, to handle infinities which cancel in 
our manipulations of it, we have introduced a temporary regulator 5 in the 
graph, which is a standard technique in the evaluation of such complicated 
integrals. It is easy to see that each of the top and bottom vertices is 
one step from uniqueness and this suggests one uses integration by parts 
on the internal vertex which includes the line with exponent (3 — fi). The 
rule for this has been given several times in previous work such as [^, 
Consequently, one obtains the difference, after one integration, of two two 
loop graphs 

(3 — [i — 5, [i — 1,/i — 1 + 5, 1,/^ — 1 — 5) — (3 — /U,/x — 1,/x— 1+5, 1,/x — 1 — 5} (5.8) 

Since the expression is multiplied by a(/i — 5) one needs the 0(5) term of 
(5.8). This is achieved by Taylor expanding each integral of (5.7) in powers 
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of 5 but since the location of 5 is common in several exponents in each term, 
expanding (5.8) gives 



— ChT(S - ii - e, 1) 



e=0 



(5.9) 



which can now be easily evaluated. Collecting terms and setting 5 to zero 
the third graph of fig. 11 is equivalent to 



7r 3 ^a(2/x - 2) 

(/x- 2)^-1; 

H 2* 

(2^-3) 

1 



a(3 - (i)a(n - 1) + 



2(M"1) 2 



1 



+ 



2fi - 3 // - 2 
1 



(2p-3)(/i-l) (/x-2)(/x-l)y (^-2)2 



2(/i - 1) 
2a 2 (l) 



+ 



+ 



(2 M -3)( M -2) 

(5.10) 



where = if)'(2(i - 1) - ^'(2 - //) - + ^'(1). This completes the 

steps required to compute ILj^- The final result is 



II5B2 



(2// - 3) 



TT b ^a(2fx - 2) 
"(/i-im/i) 

(//-!)(// -2)2 ' (/x-2) 



+ 



(/x-2) 
1 



$ + ^ 



+ 



2(/x-l)2 
27r 6 ^o 2 (2^-2) 
(^-1)6(^-2)2 



(5.11) 



For II5B1 and II6B2 a common integral lurks within each and deserves 
separate treatment. It is illustrated in fig. 12 and after the transformation 
— > one obtains the second integral of fig. 12 where we have again introduced 
a temporary regulator 5 in advance of using integration by parts on the left 
top internal vertex. This yields a set of four integrals, two of which are finite 
and proportional to the two loop graphs ChT(l,3 — fi) and ChT(l, 1) and 
two which are divergent but they arise in such a way that the 1/5 infinity 
cancels ie 



7r M a(/i-(5)a(2^-3)a 2 (l) 



a(l + 5)ChT(l,fi-l-5)- 



a(l)a((i -1 + 5) 
a(ji - 1) 



ChT(l -5,ii-l) 

(5.12) 



The finite part of (5.12) can easily be deduced by Taylor expanding each 
two loop integral so that overall the sum of contributions to the integral of 
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fig. 12 means that it is 

(2/i-3)vr 3 ^a 3 (l)a 2 (2/i-2) 



2(^-2) 



+ 



1 



66 + 



+ 



13 



2(M-1) 2 



(20 -3) 2 (2/i-3)(/x-l) (2/x-3)(/i-l 



(5.13) 



This is the hardest part of H^bi and Il6_B2 to compute. The remaining pieces 
of each can easily be reduced to two loop integrals which can be determined 
by methods we have already discussed. The final result for each is 



+ 



(2/x - 3)vr^a(2/i - 2) 

2(0- 1)5(/, -2)1^) 
1 



69 - $ 



$2 + 



2(0 - l) 2 (20 - 3) 



+ 



(20-3)(0-l) (2 M -3)(0-l) (0-1 



2(0 - 2)* + (0 - 2) 



0* 



I) 2 . 
(5.14) 



and 
1^652 



+ 



+ 



(20 - 3)7r 6 ^q(20 - 2 

(0-1)5(0-2) 

* 1 



$ f 2 3(0-1) 
2 + ~~2 ~ (20 - 3) 



G + 



(0-1) 2 



(0-2)* 



2(0 - 1) 

0* 



4(0 



I) 2 
+ 



(20-3)(0 
2 



6 



2(20 -3)(0- l) 2 ' (20-3)(0 

Calculation of 



l) 



vr 6 ^a 2 (20-2) 
(0 - 1)7(0 - 2) 



(5.15) 



There remains only one integral to evaluate. As we have already recalled 
one has to include the integral U\b of fig. 1 in order to obtain the exponent 
Ai correctly at leading order. In [Q, 11] it was determined at leading order in 
1/N. However, its 0(1/N) correction needs to be included for A2 with the 
anomalous dimensions of a, /3 and A now non-zero and we therefore define 
the 1/N expansion of the integral as 



niB = n 1B1 + 5ip + o 



1 

AT2 



(6.1) 



and UiBi = O(l). The formalism to determine IIi£2 has been discussed 
extensively in |TTJ. Basically, by using recursion relations it is possible to 
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rewrite the two bosonic integrals which occur in Hib after taking the fermion 
trace as a sum of graphs which are finite at a = fx + ^i], (3 = 1 — 77 — \ 
and A = [i — 1 + 0(1/N). As most of the integrals which then occur have 
coefficients which are 0(1/N) one can write down the contributions of these 
integrals when the fields have zero anomalous dimensions. For one integral, 
though, (a — 1, a — 2, a — 1, a — 1, £ + 2), this is not the case where £ = 3 — fi 
~ V ~ X ~ X an d we have set A = fj, — 1 + A'. We now detail its evaluation. 
Manipulation using recursion relations, [O], and using the transformations 
y and \ in the notation of j| gives 

(a — 1, a — 2, a — 1, a — 1, £ + 2) 
_ a 2 (q - l)a(g + l)(2/x - 2a - g) 
~~ 2(4a + 2£ - 3/x - 1)o(2q + f - /x - 1) 
[ 2(4a + 2g - 3fi - l)(2a + g - /x - 1) 

x . (e + i)(/u-e-2) 

x (a — 1, 2a + £ — fx — 1, 2a + £ — a — 1, 2/x — 2a — £ + 1) 
/ (3/i-4a-g + 2)(4a + g-2/i-3) \ 

v (e + i)(M-e-2) J 

x (a - 1, 2a + £ - /x - 1, 2a + £ - /i - 1, a - 1, 2/x - 2a - £ + 1}] 

(6.2) 

The coefficient of each of the two integrals can easily be expanded in powers 
of 1/N, whilst the integrals themselves need to be expanded to the same 
order. For the latter integral this is achieved by rewriting it as 

[(/x - 1 + ±r?, /x - 1 + §7/, 1 - |r/, 1 - §77, 2 - X ~ A) 
- (zx - 1 + ±77, /x - 1 + x + A - \r), 1 - ±7?, 1 - ±77, 2 - x ~ A}] 
+ (/x-l + i77,/x-l + x + A-|7?,l-i77,l-i7?,2-x-A) (6.3) 

which is an exact result where we have first of all made a conformal trans- 
formation based on the right external vertex, Q, followed by mapping the 
integral to momentum space. However, in the second term there is unique- 
ness at one of the internal vertices of integration which means it can be 
computed using fig. 5 exactly and then expanded to 0(1/N). For the first 
two terms of (6.3), since we have a difference it is easy to see that this will 
be 0(1/N). In other words we have chosen to subtract off an integral whose 
leading order value coincides with that of the integral we require, in much 
the same way as the method of subtractions is used for A divergent graphs. 
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However, here both graphs have the same structure as in fig. 8. In this lin- 
ear combination the exponents of all but one propagator coincide. Therefore 
the first two integrals of (6.3) simply become 



Xi 



Ai) 



N 



^{H-1,H- l + e,l,l,2) 



(6.4) 



e=0 



at leading order. The two loop integral can be deduced through a recursion 
relation which gives a sum of ChT(a>i,a>2) type integrals. We record 

2vr 2 ^a(l)a(2/i - 2)(2/x - 3)(p - 3) 



(H-1,»- l + e,l,l,2) 



(ix -2) 



1 



li — 2 n — 3 
3(/i-l)(/i-2)e 



- 2 



2(2/x - 3)(/i - 3) 



e + 



(6.5) 



Thus 



(a - 1, 2a + £ - // - 1, 2a + £ - li - 1, a - 1, 2li - 2a - £ + 1) 



+ 
+ 
+ 



2vr 2 ^ - 3)(2/i - 3)a(l)a(2/i - 2) 

(z""2) 
(2 M -l)( M -2) 



j _ ?7i AM - 3) 



(m-i) 

(/j-3)(2/j 2 -4// + 1 



1 - * + 

+ 



1 



1 



(jJL-l)(ji-2) 



(2/z - 3) 2(/i - 1) 
(2ll 2 - All + 1) 
(/x-3) 



3^(^ ~ 2) 
2(/i - 3) 



e + 



(6.6) 



Following a similar set of steps yields 

(a - 1, 2a + £ - // - 1, 2a + £ - //, a - 1, 2/x - 2a - £ + 1) 



tH^/j - 2)(2/x - 3)a(l)o(2/x - 2) 



x * 



2(/j - 3) 
1 1 



1 + 



771 / (2li - 1)(ll - 2) 



N 



(M-l) 



+ 



1 



(2/j - 3) 2(/j - 2) 2(/i - 1 



(2/i 2 -4 / u + l)(/i-5) (a*-1)( M -4)(2a*-5) 



2( /i -l)( M -3) 



2(//-2) 2 (/ J -3) 



(6.7) 



The remaining amount of effort in determining T1±B2 lies in simply adding 
up all the 0(l/N) terms which we believe we have done correctly due to 
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the frequent cancellation of denominator factors such as (2/i — 5), (fi — 4), 
(fi — 3) and (2/i — 3). The cancellation of the latter is reassuring as its 
potential appearance in the final answer would have indicated unwelcome 
singular behaviour in three dimensions. Overall we obtained the relatively 
simple result 



niB 



2vr 2 ^ 



(/^-i) 2 r 2 (M) 



1-^ 



3/x(2/x - 3) 



G + 



1 



(6.8) 



where we record that we have used the results 



(fi - 3, n - 1, fi - 1, n - 1, 5 - fi) 



a(5 — fi) 



(fi - 3) 2 C7iT(l, 1) + 



7r^(/i 3 - 10/i 2 + 31/i - 31)a(l) 



a(l)a 2 (2) 
(/i - 2, /j - 1, n - 1, /x - 2, 5 - fi) 



(fi - 2)3 a (3 - fi) 



a(5- fi)(fi 
a»(2) 



C7iT(l,l) 



^(/i - l)o(l) 
(^-2)3a(3-M) 



(6.9) 



(fi - 2, fi - 1, fi - 2, fi - 1,5 - ft) 



a(5 — /i) 



C7iT(l,l) + 



tt 2 ^(2/j 2 - 12/J + 19)o(2) 
(A* - 2)3a(3 - /x) 



which were incorrectly given in [JO]], but were not required for that work. 
We close our discussion of the evaluation of our graphs by noting that the 
correct expressions for the integral F(a, (5) defined in is 



F(a,(3) 



+ 



ir 2 ^a(fi - l)a(2fi - 2) 
(a + (3 + fi - S)(2 - a - 0) 



(v-iy 

which is valid for all a and j3. 



a(a)a(2 - a) + a(f3)a(2 - 0)] 
ChT(a,f3) (6.10) 



7 Discussion. 



The previous three sections have been devoted to the evaluation of the inte- 
grals which appear in the formal master equation (3.7) and (3.8). It is now 
a straightforward matter of substituting for the various expressions in each 
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equation and evaluating the 0(1/N 2 ) correction to the determinant. As an 
intermediate step we note, 



n 2B 



7T 



+ 



+ 



- i) 2 r 4 (M) 

16 

10 

+ 



8(2// - 3) 



$ + ^ 2 ) 



4(2/i-3)(/x-3) 
(/i- 1)(m-2) 2 
2 16 



+ 



12(2/i - 1)9 
G""2) 

16 : 



(/x-1) 2 (/x-2) 



(/i - 1) (/x - 2)2 ^ - 2)2 m 
and also record that 

/xr 2 (/x)r7 2 



^2 



2tt 2 m(^- 1) 



(M-l) 



+ 2 + (2/i -!)*(//) 



(7.1) 



(7.2) 



With these expressions together with the expansions of p(P), q((3), r(a — 1) 
and s(a — 1) to the next to leading order values, we find from the vanishing 
of the determinant of the matrix defined by (3.7) and (3.8) at 0(1/ N 2 ) 



A, 



+ 



+ 



(2/i - 3)/x 



Qi - 2)2 



+ 



1 



(M-2) 
.2 3 



($ + ^ 2 



3^e 



7 



9 - 2/x + 



2(/x - 2 
6 



/x- 2 



1 



+ 



■ - 2// z 

+ 2/x 2 
1 



2 

5/i 



1 

2/i 

■3 + 



4/i 



-i; 
i 

v 



2(/i-l) (/x-l)2 4(/i-2) 2(/x-2)2 



(7.3) 



which is an arbitrary dimensional expression for the 0(1 /N 2 ) corrections to 
the /3-function of (2.2). It is worth recording that an independent check on 
the correctness of (7.3) is that it ought to agree with the expansion of the 
critical /3-function slope computed explicitly from the three loop result of 
(2.2). We have checked that this is indeed the case. Further, we can deduce 
the value of the exponent in three dimensions as 



A 



16 32(27tt 2 + 632) 



37r 2 AT 



+ 



27vr 4 iV2 



(7.4) 



As two independent exponents v = 1/(2A) and rj or rj + x are now known 
at 0(1/ N 2 ) this implies that the remaining thermodynamic exponents of 



21 



the model can be deduced through the hyperscaling relations which were 
recently checked at leading order in (l3j . 

With (7.4) we can now gain an improved estimate of the exponent v 
in three dimensions and compare with recent lattice simulations where the 
same exponent is calculated for the case = 8 in our notation, [14|. Thus 



1 + 



32 64(27vr 2 + 584) 



(7.5) 



3vr 2 A 27tt 4 7V 2 
and employing a Pade-Borel technique widely used in improving estimates 



of exponents, 21, 29, pCl, we have 



N 



dt e 



-Nt 



m t 32(27vr 2 + 608)t 2 
3^ ' 



+ 



81tH 



0.98 (7.6) 



N = 8 



Recent simulations, 14], give v = 0.98(7) and so (7.6) is in excellent agree- 
ment with that Monte Carlo result. Also the exponent 2 — r] — x, in our 



notation, has been calculated numerically as 1.26(3) in [14] and we record 
that from (2.6)-(2.9), at N = 8, we find 2 — r] — x = 1-25 again in good 
agreement. 

We conclude by making several remarks. First, the Gross Neveu model 
has now been solved at 0(1/N 2 ) The techniques we have had to employ are 
different from those used to perform the analogous calculation in the 0(N) a 
model, due to the appearance of several four loop graphs. More importantly, 
though, we have laid a substantial amount of the groundwork for performing 
the same calculation for QED. Whilst this is a more complicated theory the 
basic techniques to treat the integrals have been developed here. 



Acknowledgement. The author thanks Leo Karkkainen and Pierre La- 
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Note added. Whilst in the final stages of this work we received a preprint, 
|j3l[ , where A2 is stated and we record that it and (7.3) are in agreement. 
We believe the method of |H]] is different from the one given here. 
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Figure Captions. 

Fig. 1. Skeleton Dyson equations to determine Ai. 

Fig. 2. Higher order graph for if) consistency equation. 

Fig. 3. Divergent higher order corrections to a consistency equation. 

Fig. 4. Finite higher order graphs for a consistency equation. 

Fig. 5. Uniqueness rule for a bosonic vertex. 

Fig. 6. Uniqueness rule for a fermionic vertex. 

Fig. 7. Subtractions for ~R-2B\ and ~Rzb- 

Fig. 8. Definition of {ai, ^2, 03, 04, 05} and (ai, 0:2, 03, 0:4, 0:5). 

Fig. 9. Intermediate integrals in the evaluation of ILlb- 

Fig. 10. Il6£i after several integrations. 

Fig. 11. Intermediate integrals in the evaluation of ~R<oB2- 

Fig. 12. Common integrals in the evaluation of IT5B1 and II6B2- 
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